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ABSTRACT 

Context. Models of primordial and hyper-metal-poor stars that have masses similar to the Sun are known to experience an ingestion of 
protons into the hot core during the core helium flash phase at the end of their red giant branch evolution. This produces a concurrent 
secondary flash powered by hydrogen burning that gives rise to further nucleosynthesis in the core. 

Aims. We aim to model the nucleosynthesis occurring during the proton ingestion event to ascertain if any significant neutron-capture 
nucleosynthesis occurs. 

Methods. We perform post-process nucleosynthesis calculations on a one-dimensional stellar evolution calculation of a star with mass 
1 M and a metallicity of [Fe/H] = -6.5 that suffers a proton ingestion episode. Our network includes 320 nuclear species and 2,366 
reactions and treats mixing and burning simultaneously. 

Results. We find that the mixing and burning of protons into the hot convective core leads to the production of 13 C, which then burns 
via the l3 C(a,n) l6 reaction releasing a large number of free neutrons. During the first two years of neutron production the neutron 
poison 14 N abundance is low, allowing the prodigious production of heavy elements such as strontium, barium, and lead via slow 
neutron captures (the s process). These nucleosynthetic products are later mixed to the stellar surface and ejected via stellar winds. 
We compare our results with observations of the hyper-metal-poor halo star HE 1327-2326, which shows a strong Sr overabundance. 
Conclusions. Our model provides the possibility of self-consistently explaining the Sr overabundance in HE 1327-2326 together with 
its C, N, and O overabundances (all within a factor of ~ 4) if the material were heavily diluted, for example, via mass transfer in a 
wide binary system. The model produces at least 18 times too much Ba than observed, but this may be within the large modelling 
uncertainties. In this scenario, binary systems of low mass must have formed in the early Universe. If true then this puts constraints 
on the primordial initial mass function. 

Key words, nuclear reactions, nucleosynthesis, abundances - stars: evolution - stars: individual (HE 1327-2326) - stars: interiors - 
stars: Population II - stars: Population III 



1. Introduction 

Stars with a similar mass to the Sun eventually expand and be- 
come red giants and at the end of their red giant phase undergo 
a thermonuclear runaway in the core known as the core helium 
flash. The ignition of helium in this flash event occurs at the 
hottest place in the core, which is located off centre due to cool- 
ing via neutrino energy losses deep in the stellar interior. The 
amount of energy released by the He flash is so large that it can 
not be carried radiatively, and thus convective motions are trig- 
gered, forming the so-called He convective zone (HeCZ). A typ- 
ical HeCZ extends from the location of the He-burning region 
to close to the base of the H-burning shell. The convection does 
not normally break through the H-burning shell because the ac- 
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tive H-burning shell provides an entropy barrier against mixing 
dlben Jr.lll976h . 

However, in solar-mass red giant stars of primordial (Z = 0) 
or hyper- metal-poor ([Fe/H] < -5.0) composition the H-He in- 
terface is not as impenetrable as in more metal-rich red giant 
stars. The reason for this is twofold. First, the core helium flash 
starts much further off centre in these low-metallicity models (at 
~ 0.27 M in the model we present here, Fig. [TJ than at so- 
lar metallicity and is thus closer to the H shell. Second, the en- 
tropy barrier at the H shell is much weaker in stars of very low 
metallicity because the H-burning shell alm ost switches off at 
this stage of evolution (Fujimotoetal. 1990). Thus, in primor- 
dial solar-mass red giants the He convective region can pene- 
trate into the overlying H-rich region and protons can be car- 
ried down into the hot He-burning zone. If enough protons are 
ingested a secondary (H-burning) "flash" will occur within the 
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Fig. 1. Internal evolution of our solar-mass model star with a 
metallicity of [Fe/H] = -6.5 during the core He flash. Proton 
ingestion starting at t — (6.3 Gyr into the stellar evolution) due 
to the expansion of the HeCZ leads to the production of 13 C and 
the consequent neutron superburst via 13 C(a,n) 16 0. Convective 
regions are grey, HeCZ: He convective zone, HCZ: H convective 
zone. The blue solid line and the black dashed line represent the 
location in mass where the abundance of H is 0.0001 and 0.73, 
respectively. The insert is a zoomed-out view to help localize the 
region of interest within the star as a whole. It can be seen that 
some of the processed material starts to be dredged up by the 
encroaching convective envelope at around t = 1000 yr. 



He convective region while the He-burning flash is still ongo- 
ing. This can result in the He-convective region splitting into 
two convective regions at the location where the H burning re- 
leases most of its energy. This sequence of events is well estab- 
lished in ID stellar models and is referred to variously as "dual 
core flash" (DCF), "proton ingestion ep i sode" or "helium flash 
induced mixi ng" (L Fuiimoto et al. 1990; Hollowell et al. 1990; 
Schlattl et all 120011; IPicardi et alj |20o4 Icampbell & Latta nzio 



2008t ISuda & Fuiimotoll2010l) . Importantly, the material pro- 



cessed by the DCF is later mixed to the stellar surface by a (sin- 
gle) dredge-up event. Thus the nucleosynthetic products of the 
DCF are available to be shed into the surrounding interstellar 
medium by stellar winds or transferred to a binary companion. 
The main chemical products are C and N, and this qualitatively 
agrees with the observations of C-rich extremely metal-poor 
halo s tars in the Galactic Halo (ISchlattl et al.ll200 it IPicardi et all 
l200l . 

Here we explore in detail the consequences of a proton in- 
gestion episode during the core He flash on the nucleosynthesis 
and surface abundances of a star of 1 M Q and [Fe/H] = -6.5. 



2. Method 

We have used one of t he stel lar structure models computed in 
Campbell & Lattanzio! (2008) (Paper I). This model was cal- 
culated using the Schwarzschild criterion to define the border 
of convective regions and a diffusive mixing treatment of the 
m ovement of material inside convective regions similar to that 
of iMevnet et alj d2004l) . The stellar structure code follows six 
chemical species: 'H, 3 He, 4 He, 12 C, 14 N, and le O. The rest of 
the isotopes involved in the CNO cycles are included by assum- 
ing that they are always present in their equilibrium abundances. 

We calculate the detailed nucleosynthesis using a post- 
processing code that activates nuclear reactions in the star on the 



basis of the information on the temperature, density, and con- 
vective velocities provided by the stellar structure model (see 
Paper I and references therein for details). The equations that 
describe the changes in the abundances are solved using an im- 
plicit method so that a large matrix of n 2 is solved, where n is the 
number of species in the network. Nuclear burning and convec- 
tive mixing are both included in the equations of the abundance 
changes and thus solved simultaneously. Our present calcula- 
tion was performed using a nuclear network of n — 320 species, 
from protons and neutrons up to lead and bismuth, and 2,336 re- 
action s with rates from the JINA reaclib database ( Cvbur t et al.1 
120101) . The initial composition for the elements up to Zn was 
taken as a mix of standard Big Bang material and supernova 
calculations as described in Paper I. The abundances of the el- 
ements heavier than Zn were taken equal to zero, however we 
note that this choice is mostly irrelevant to our results because 
the neutron-capture nucleosynthesis overwrites the memory of 
the initial abundances. 

Due to our computationally expensive scheme we only in- 
cluded some of the s-process branching points that may be acti- 
vated under conditions of a high neutron flux. This shortcoming 
should not have a strong impact on the overall heavy-element 
distribution reported here because the overall neutron flux is only 
very marginally affected by the details of the s-process path. 

Due to computational problems we have not yet proceeded 
much past the end of the neutron superburst. However, this did 
not prevent us from accurately estimating the surface composi- 
tion of the neutron-capture elements for the model (Sec. 3.2). 
This is because there are no further mixing episodes following 
the dredge-up event that mixes the products of the core He flash 
into the stellar envelope. 



3. Results 

3.1. The "Neutron Superburst" 

We find that our detailed DCF model experiences a phase of very 
high neutron flux in the He-rich convective region, which we call 
a neutron "superburst". The neutron flux for the superburst is pri- 
marily supplied by 13 C(ar,n) O reactions occurring at the base 
of the HeCZ. The 13 C for this reaction is produced as a conse- 
quence of the proton ingestion - it is a product of partial CN- 
cycling via 12 C(p,y) 13 N(yS + ) 13 C, where the 12 C is supplied by 
the on-going triple- g reactions of t he He-fla sh (Campbell 2007; 
lLugaro etal!l2009l) . We note that iFuiimoto et all (ll990h spec- 
ulated that neutrons and light s-process elements may be pro- 
duced during this event. 

The superburst occurs during the earliest phases of the inges- 
tion (Fig. [TJ. This is because the H profile is very extended (Fig. 
[2) due to the nature of the previous H-shell burning. In a star with 
a negligible amount of CNO catalysts the p-p chain reactions 
dominate over the CNO cycles in terms of the energy produc- 
tion. Since the p-p chains have a low temperature dependence, 
the H-burning shell is active over an extended region, leading to 
a broadened H-abundance profile. The presence of this long H 
"tail" leads to variations in the amount of protons ingested into 
the expanding HeCZ. In the early ingestion phases, the mixing 
of small amounts of protons leads mostly to production of 13 C 
nuclei and the consequent neutron superburst (t = to ~ 1 .2 yr 
in Fig. 0. After about 1.7 yr, the HeCZ has expanded through 
the entire H abundance profile up to the base of the unburnt 
envelope. At this point the abundance of protons in the HeCZ 
becomes so high that the abundance of the neutron poison 14 N 
exceeds that of 13 C owing to complete CN-cycling. Under these 
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Fig. 2. The H profile around the location of the proton ingestion, 
at various times during the neutron superburst (f = corresponds 
to the same time in Fig.Q}. The peak neutron flux occurs between 
t ~ and t ~ 1.2 yr. 
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Fig. 3. Evolution of the neutron density (lower panel, in loga- 
rithm) and of the abundances of selected heavy element species 
(upper panel, in logarithm) as a function of time during the main 
phase of the neutron superburst (t — corresponds to the same 
time as in Figs. [1] and [2] and to the onset of the neutron flux. 
Abundances are sampled at 0.29 M , the location of the maxi- 
mum neutron density. All abundances, except those of protons 
and neutrons, are homogeneous within the HeCZ, and hence 
their values in the plot are representative of the abundances in 
the whole region. 



conditions the 14 N(n,p) 14 C reactions dominate and the neutron 
superburst is quenched. 

The neutron density during the superburst peaks at ~ 10 15 
n/cm 3 , remains above 10 14 n/cm 3 for about 0.2 years (Fig. [3]), 
and the time integrated neutron flux is 287 mbarri -1 . At the loca- 
tion of the peak neutron density, at 0.29 M Q in mass, and for neu- 
tron densities > 10 12 n/cm 3 , the temperature decreases in time 
from 212 to 160 MK and the gas density from 15,200 to 2,900 
g/cm 3 . Free neutrons are available for neutron captures that re- 
sult in efficient production of the elements heavier than iron and 
up to lead (Fig. [3}. The large total neutron flux in combination 
with the long time scale for the release of the neutrons result in 
elemental abundances quite typical of the s process at low metal- 
licity. We find production of nuclei at the three s-process peaks 
corresponding to the magic number of neutrons 50, 82, and 126, 
as represented by the isotopes of 88 Sr, 138 Ba, and 208 Pb in Fig. [3] 

At the start of the neutron burst the abundance of 56 Fe is de- 
pleted as the neutron capture chain converts the iron into 88 Sr. 
After the neutron exposure reaches the value of =2 10 mbarrT 1 , 
the abundance of 88 Sr starts to decrease and that of 138 Ba to in- 
crease. This is caused by neutrons overcoming the first bottle- 
neck at the magic neutron number N = 50. Slightly later on, as 
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Table 1. Selected elemental [X/Fe] = log l0 (X/Fe) sta r • 
except for Li given as log e(Li) = log (N Li /N H ) + 12). 
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the neutron exposure reaches ^ 20 mbarn -1 the same applies to 
the abundance of I38 Ba as the second bottleneck at TV = 82 is 
overcome and 2()8 Pb starts being produced. Once the neutron ex- 
posure reaches above ^100 mbarn -1 the absolute abundances 
increase until the end of the neutron flux. By the end of the neu- 
tron flux the abundance of 56 Fe has grown by a factor of 232 
due to neutron captures on the lighter elements starting from the 
abundant 12 C. 



3.2. Surface composition and comparison to the halo star 
HE 1327-2326 

Dredge-up of the material previously included in the HeCZ oc- 
curs about 10 4 years later, where we find the maximum post- 
flash penetration of the convective envelope. In Table Q] (Column 
2) we show the resultant surface abundances of the model after 
this dredge-up event. Owing to numerical problems, our current 
calculations only reached to a time just after the splitting of the 
convective zone. Thus, in Table[T]the surface abundances of the 
elements up to Al have been taken from Paper because they 
are affected by further proton captures after the convective zone 
splits, which were not included in the present calculations. The 
abundances of the elements heavier than Al are evaluated by 
considering that ^0.1 M of the material previously involved 
in the neutron superburst is dredged-up to the stellar surface 
and diluted with an envelope mass of 0.44 M Q , corresponding 
to a dilution factor f dredge-up - 0.23. The surface abundances 
are then calculated for each nuclear species as mass fractions: 

X surface ~ [^initial + X[]eCZ X f dredge— up\ I [ 1 fdredge-up\j where 

^initial is the initial envelope value and Xfj e cz is the value in the 
HeCZ. 



1 This is self-consistent because the nuclear species included in our 
present extended network have so much smaller abundances than the 
lighter elements that their effect on the overall neutron flux is insignifi- 
cant. 
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It can be seen in Table Q] that neutron-capture elements ap- 
pear at the surface. Their abundances are highly enriched, in the 
case of Pb reaching up to the level of the absolute solar abun- 
dance. The high neutron density of the neutron flux contributes 
to a distribution that is weighed towards neutron-rich nuclei. For 
example, 87 Rb is produced resulting in Rb/Sr =s solar. Fe is also 
increased by 2.2 dex. As reported in Paper I, the model also pre- 
dicts a large overproduction of C, N, O, and F and milder en- 
hancements of Na and Mg. 

In Column 3 of Table Q] we include the composition result- 
ing from further diluting the surface composition with material 
of initial composition by a factor /oil — 3 x 1CT 4 . This specific 
dilution factor was used to match the C abundance derived from 
the most metal-poor star discovered to date, the subgiant star 
HE 1327-2326, with [F e/H] = -5.96 (Column 4 of Tab leffl with 
typical error bars +0.25: lAoki et al.l2006UFrebel et all 2008). By 
applying this further dilution factor our model self-consistently 
reproduces the observed C, N, and O abundances within a factor 
of 4 for N and 3 for O (observational error bars are ~ a factor 
of 2) with Fe increasing by only 20%. However, it produces 59 
times more Sr and at least 260 times more Ba than observed. 
We performed a preliminary post-processing calculation using 
an initial [Fe/H] = -5.76 (numbers in brackets and italics in 
Column 3 of TableQ}, more similar to HE 1327-2326. We would 
not expect such relatively small change in the Fe abundance to 
have a major impact on the overall stellar structure and the nu- 
cleosynthesis of the light elements, but we discovered that it has 
a strong impact on the nett yields of the neutron-capture process, 
while the relative abundance distribution is not significantly al- 
tered. The results from this calculation show a better match to 
the observed Sr abundance, within a factor of 4, while still pro- 
ducing too much Ba by a factor of at least 18. In this case, Fe is 
unchanged. 

The composition of HE 1327-2326 could have arisen from 
binary system mass-transfer via wind accretion or Roche-lobe 
overflow from a star such as that modelled here (which would 
now be a white dwarf). Subsequent dilution of the accreted ma- 
terial in the envelope of HE 1 327- 2326 via convection or thermo - 
haline mixing would be expected (Stancliffe & Glebbeek 2008). 
Since radial velocity variations have not so far been detected for 
HE 1327-2326, the Roche-lobe overflow mass transfer scenario 
is not favoured. Wind accretion in a very wide binary remains 
a possibility and would be more in line with the large dilution 
factor that we have estimated. The amount of mass that a binary 
companion of mass 0.8 M would need to accrete from our hy- 
pothetical primary star to match the observed C abundance is 
~ 2.4 x 10 -4 M . Here we have assumed the extreme case where 
the accreted material mixes throughout the entire star. Using this 
accretion mass together with an estimated total mass loss from 
the primary star of 0.4 M (leaving a 0.6 M H WD) the expecte d 
period of the system is 769 yr, using Eq. 12 of Su da et all {2004), 
which includes Bondi-Hoyle accretion. This is consistent with 
the current non-detection of radial velocity variations. 

4. Discussion 

We compared our theoretical surface abundances with the ob- 
served abundances of HE 1327-2326. Due to the high degree 
of surface pollution in our model only a small amount of mass 
should have been accreted by HE 1327-2326 to match its com- 
position, indicating a wide binary configuration. The model self- 
consistently reproduced C, N, O, and Sr within a factor of 4 (if 
[Fe/H] = -5.76), but overproduced Ba. The model produces 
more Sr than Ba before 0.008 yr (Fig. [3]), when the neutron 



exposure is less than ^10 mbarri -1 . At the first peak of 88 Sr 
we find [Ba/Sr] < 0.19, which would be consistent with the 
observed [Ba/Sr] < 0.23. However, the absolute Sr abundance 
in the HeCZ at this point is 120 times smaller than the final 
abundance dredged-up (for the [Fe/H] = -5.76 model), while 
a factor of =* 17 smaller could still provide a match to Sr in 
HE 1327-2326 within a factor of 4. To match the composition 
of HE 1327-2326, in this case also including Ba, we would need 
C, N, and O - 7 times smaller and fan — 7 times larger than 
reported in TableQ]with an expected binary period of 179 yr, as- 
suming the same amount of dredge-up. These possibilities may 
be within the model uncertainties and need to be carefully in- 
vestigated as the DCF scenario appears to have the potential to 
self-consistently reproduce m ost of the strange overabundances 
in hyper-metal-poor stars (e.g.lSchlattl et aljfeOQlt iPicardi et al.l 
120041: ISuda & Fuiimotol l2010l ; Paper I; Table ffll. The DCF 
would be best studied in the framework of multidimensional 
models (lDeupreell996r;lDearborn et alj2006tlMocak et al.l2010t 
iHerwig et alkoiOl) and effects such as convective overshoot, ex- 
tra mixing, and rotationally induced mixing should be also inves- 
tigated in this context. 

Our model underproduces Na, Mg, Al, Ca, Ti, and Ni. This is 
suggestive that these elements came with F e from an early super - 
nova that polluted the protostellar cloud dJoggerst et al.l 12010). 
even though Na is also underproduced in these models. The un- 
certainties of the Ne and Na proton-capture reactions must be 
investigated. Lithium is another mismatch since the DCF low- 
ers the surface Li abundance by only a factor of two while 
HE 1327-2326 is heavily depleted in Li. This problem needs to 
be addressed in terms of mixing processes on the secondary star 
HE 1327-2326. 

The peculiar chemical signatures of the DCF event may be 
able to shed light on star formation and the initial mass function 
in the early Universe, since the implication of our scenario is that 
binary systems of roughly solar-mass stars should have formed 
from hyper-metal-poor gas. We suggest that our results should be 
taken into account in the study of the composition of extremely 
metal-poor stars in our halo and nearby galaxies. 
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